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Summary. — Statistical mechanics is one of the most comprehensive theories in 
physics. From a boiling pot of water to the complex dynamics of quantum many- 
body systems it provides a successful connection between the microscopic dynamics 
of atoms and molecules and the macroscopic properties of matter. However, statis¬ 
tical mechanics only describes the thermal equilibrium situation of a system, and 
there is no general framework to describe how equilibrium is reached or under which 
circumstances it can be reached at all. This problem is particularly challenging in 
quantum mechanics, where unitarity appears to render the very concept of thermal- 
ization counterintuitive. With the rapid experimental progress in the control and 
probing of ultracold quantum gases this question has become within reach of de¬ 
tailed experimental investigations. In these notes we present a series of experiments 
with ultracold one-dimensional Bose gases, which provide novel insights into this 
fundamental question. 
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1. Motivation 

Does an isolated quantum many-body system that is prepared in a non-thermal initial 
state relax to thermal equilibrium? As we know from our everyday experience, many 
physical systems can very successfully be described by a thermal state. On the other 
hand, the time-reversal symmetry that results from the unitarity of quantum mechanics 
seems to make the relaxation to thermal states impossible in an isolated system [1], 
This seemingly simple question thus addresses the fundamental relation between the 
macroscopic description of statistical mechanics and the microscopic quantum world. It 
has been highly contested since the 1920s [2] and important theoretical advances have 
been achieved over the years [3, 4, 1, 5]. Variations of this question play important 
roles in such diverse helds as cosmology, high-energy physics and condensed matter [6, 
7, 8, 9, 10, 11]. However, it has only been through the recent experimental progress 
in manipulation of ultracold quantum gases that this question has become within reach 
of detailed experimental investigations [12]. In the following we will present a series 
of such experiments, which we performed using ultracold one-dimensional Bose gases. 
The versatility of these gases allowed us to realize several textbook-like non-equilibrium 
phenomena, which provide important insights into the dynamics of quantum many-body 
systems. 

2. One-dimensional Bose gases 

Over the last years, one-dimensional (ID) Bose gases have proven to be a versatile 
testbed for the study of quantum many-body physics in and out of equilibrium. The 
great interest in these systems stems from several key properties. From the theorist’s 
perspective ID Bose gases offer a rich variety of interesting many-body effects, while still 
being tractable with reasonable effort [13, 14]. On the experimental side their realization 
using cold atomic gases offers precise control over many system parameters, as well as 
highly-effective means to probe their dynamics [15]. In this first chapter, we will briefly 
outline important aspects of the description of ID Bose gases. For more detailed accounts 
we refer the reader to Refs. [13, 16, 17]. 

The experimental realization of a ID Bose gas follows the familiar procedure based 
on laser and evaporative cooling that is also used for the production of Bose-Einstein 
condensates from three-dimensional (3D) Bose gases [18, 19]. However, creating an ef- 
fectivly ID system in a 3D world requires extremely asymetric traps with a very tight 
confinement in all but one spatial directions. The general aim of this tight confinement is 
to raise the energy splitting of ground and first excited state in the two tightly-conhned 
directions, such that all relevant energy scales of the trapped gas lie below it. For a 
harmonic trap this means that the temperature T and the chemical potential fj, fulfill 
ksT,^ ^ Huj±, with ks denoting Boltzmann’s and h the reduced Planck constant. This 
realizes a situation where the dynamics along the radial directions can be integrated out 
leaving the dynamics along the weakly conhned axial direction described by an effective 
ID model. Contact interactions in this ID model can be parametrized by an effective 
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scattering potential with the interaction strength [20] 

(1) g = 2hasU!±. 

Here, o* is the s-wave scattering length of the gas. Note that this description assumes 
that microscopic scattering processes still have a 3D character, which is the case as 
long as the s-wave scattering length is small compared to the ground state width of the 
tight radial confinement, i.e. Us yJhJmwZ, with m denoting the mass of the atoms. 
Interesting effects like confinement-induced resonances can occur when this assumption 
is no longer valid [20, 21]. 

Such highly-anisotropic trap conhgurations can be created in strongly-focussed optical 
dipole traps [22, 23, 24], optical lattices [25, 26, 27, 15] or in magnetic micro traps [28, 29]. 
In our experiments, we rely on the latter because micro traps, as we will see below, 
allow for a particularly precise and convenient preparation of non-equilibrium states. 
Typical trap frequencies in our setup are oj±_ = 27r • 2 kHz in the tightly-confining radial 
directions and Wax = 27r • 10 Hz in the weakly-confining axial direction. The ID Bose gas 
is created in this trap by evaporative cooling of an elongated 3D thermal cloud through 
the condensation crossover and then further into the ID regime. 

While the preparation of an ultracold ID Bose gas is similar to the one of an ultracold 
3D Bose gas, significantly different physics arise once the gas enters the ID regime. The 
Mermin-Wagner theorem [30] tells us that no off-diagonal long-range order can emerge 
due to the enhanced role of fluctuations in ID. Consequently, there is no macroscopic 
occupation of the lowest momentum mode even at T = 0. Thus no true Bose-Einstein 
condensation is possible. Instead a large number of distinct degenerate regimes emerges 
[31, 32], which might or might not share some of the familiar features of a Bose-Einstein 
condensate. 

In the homogeneous limit the system is described by the Lieb-Lininger Hamilto¬ 
nian [33] 


( 2 ) 



dz 


d^{z) 

dz dz 




where the 'I'(z) denote bosonic field operators. The Lieb-Lininger Hamiltonian is a prime 
example of a so called integrable model [33, 34, 35, 36]. Such models are characterized 
by a large number of conserved quantities and have historically been an important topic 
in mathematical physics. Experiments with ID Bose gases can thus provide a link be¬ 
tween the corresponding deep mathematical insights and physical reality. Most notably, 
the conserved quantities have a profound influence on the non-equilibrium dynamics of 
these systems, which makes them particularly interesting for the study of relaxation and 
thermalization processes [37, 38]. 
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The interaction strength in Eq. (2) can be parameterized by the Lieb-Lininger pa¬ 
rameter 7 = mg/h'^nid- Notably the interaction strength increases for decreasing par¬ 
ticle densities nid- For 7 1 the gas is in the strongly-interacting Tonks-Girardeau 

regime [25, 39]. All experiments presented in these notes are performed with 7 ^ 1 , 
where the gas is a weakly interacting quasi-condensate. In this regime density fluctua¬ 
tions are suppressed and the density distribution is similar to the Thomas-Fermi profile 
of a BEG. However, the phase fluctuates strongly along the length of the system. 

The suppression of density fluctuations allows us to employ a generalized version of 
the well-known Bogolibov expansion even though there is no macroscopically occupied 
mode [40]. To that end, we express the field operators in terms of density and phase 
operators 

(3) #( 2 ) = 


which satisfy the bosonic commutation relation 

(4) [h{z), 9{z')] = i5{z - z'). 

Inserting this definition into the Hamiltonian in Eq. (2) leads to a quadratic model 
describing the low-energy limit of the system. The result is known as the Luttinger 
liquid Hamiltonian 


( 5 ) 



dz 


K / d9{z) 

TT \ idz 


+ ln{zf 


k 


The parameters in this Hamiltonian are the speed of sound c = ^/gnidjm and the 
Luttinger parameter K = nidihn)^/Agm. The corresponding eigenmodes are non¬ 
interacting phonons with momentum k, linear dispersion relation tok = ck and energies 
hujk- The creation and annihilation operators dk and aj, define the phonon occupation 
number fik = d\dk- They are directly related to the Fourier components of density and 
phase via 

(6) hfe ~ [ak{t) + 9k ~ [dk{t) - ■ 

One therefore also speaks of the phase and density quadrature of a phonon. Finally, we 
note that, besides cold atoms, the Luttinger liquid Hamiltonian also plays an important 
role in both bosonic and fermionic condensed matter systems [41, 42, 43, 44]. 


3. Creating a non-equilibrium state 

As we have already noted above, the main tool in all experiments that are presented 
in these notes is a magnetic micro trap. This micro trap is realized using an atom 
chip [29], a collection of current-carrying gold wires, which are micro-fabricated onto a 
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silicon substrate. Apart from the possibility to create traps with the necessary aspect 
ratio to reach the ID regime (magnetic field gradients scale as 1/r^ with the distance 
r to the current carrying structure and micro traps allow the positioning of the atoms 
at very small distances r ~ 100/rm), the atom chip also allows for a precise dynamical 
control over the trap parameters. For example, the initial harmonic trap can transversely 
be transformed into a double well potential. This is realized by radio-frequency (RF) 
dressing of the magnetic sub-states of the atoms [45]. The RF fields are applied through 
additional wires on the chip, which due to their proximity to the atoms allows for very 
high RF field amplitudes and a precise control over the field polarization. 

We use this technique to coherently split a single ID Bose gas into two halves, thereby 
creating a non-equilibrium state [46, 47]. The process of splitting is performed fast 
compared to the axial dynamics in the system so that tgpiit < ^h/c = h/fi. Here = 
U/mc is the healing length, c = the speed of sound and /x the chemical potential. 

The fast splitting assures that no correlations can build up along the axial direction such 
that the splitting happens independently at each point in the gas. The process can be 
intuitively pictured as a local beam splitter where each atom is independently distributed 
into the left or right half of the new system. The corresponding probability distribution 
for the local number of particles N on each side is therefore binomial 

(7) PiNi,Nr) = 

with Pi = 1/2 for a balanced splitting process. The resulting fluctuations in one half of the 
system are thus given by Var[A?pr] = NPi (1 —pi), which translates into (|AA^|^) = A?/4 
for AN = {Ni — Nr)/2 in the balanced case. Figure 1 illustrates this process. 

Once we can speak of two spatially separated systems we can perform a variable 
transformation to anti-symmetric and symmetric degrees of freedom, which will help us 
to better describe the quantum state after the splitting. In the following these will also 
be referred to as relative and common degrees of freedom. Starting from the density and 
phase fluctuations in the left and right halves (denoted by hi^r{z) and 6i^r{z), respectively) 
we find 


(8) ${z) = ^z) - 0i{z) , 

for the phase, and 

X flriz) — hl(z) ^ / X ^ / X 

(9) v{z) = - - - , Vcomiz) = nr{z) + ni{z) 

for the density. The usefulness of this approach becomes clear as we return to the shot 
noise, which now only enters in the relative number fluctuations 


( 10 ) 


{v{z)v{z')) 
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Fig. 1. - Local number fluctuations, (a) Schematic representation of number and phase 
fluctuations in a ID Bose gas. The boxes indicate a course graining on the length scale of the 
healing length, (b) The splitting distributes the atoms on each of these grid points binomially 
between the two wells. This results in number fluctuations with a variance of N/A (see text) in 
each gas. These fluctuations correspond to an energy which is added to the relative degrees of 
freedom of the system during the splitting. Figure adapted from [16]. 


Here, nid denotes the mean density in a single gas after splitting, which results in the 
additional factor of 2 as compared to the binomial fluctuations that were introduced 
above. Transforming these fluctuations into momentum space gives 

( 11 ) {vkVk') = -^5k-k'- 

From the commutation relation in Eq. (4), we see that the corresponding shotnoise 
introduced to the phase quadrature of the relative modes goes with 1 /nid and is therefore 
negligible. 

Returning to the Luttinger Hamiltonian (Eq. (5)) we can identify the amount of 
energy that is introduced into each individual phononic mode during the splitting process 
as 5 nid/ 2 , which is typically significantly smaller than the thermal energy of the initial 
gas. Moreover, as we have just shown this energy is only stored in the density quadrature 
of the relative degrees of freedom, while it should be equipartitioned between phase and 
density quadrature in thermal equilibrium. 

The situation is different for the common degrees of freedom, which inherit all thermal 
excitations that were present in the initial gas before the splitting. The state created 
by splitting is thus also out of equilibrium in this respect, as the common degrees of 
freedom contain a lot of thermal energy, while the relative degrees of freedom only contain 
quantum shotnoise. 

In experiment, the equilibrium situation can be realized by the transforming the 
harmonic trap into a double well while the gas is still thermal. Further independent 
evaporative cooling in both wells then results in two degenerate gases with no knowledge 
of each other, which corresponds exactly to thermal equilibrium. The experiment thus 
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enables the unique possibility to contrast non-equilibrium and thermal states in identical 
settings. 

4. Probing the quantum state 

Information about the system and its dynamics after the splitting is extracted using 
standard absorption imaging [48] after releasing the system from the trap. If only a 
single gas is present it simply expands in time-of-flight (TOF), while a pair of condensates 
expands, overlaps and forms a matter-wave interference pattern [17]. The resulting cloud 
is subsequently illuminated by a resonant laser beam, casting a shadow that is imaged 
onto a CCD camera. This method is destructive, therefore many identical realizations 
are necessary to probe a time evolution. It is important to note that the tight transversal 
confinement of the ID gases leads to a very rapid radial expansion, which results in an 
immediate dilution of the system. Therefore, interaction effects in the expansion are 
negligible and the absorption images enable comprehensive insights into the properties 
of the initial trapped system. 

A schematic overview of imaging probes employed in our experiment is shown in 
figure 2. In the following we will give a short overview of the insights into the dynamics 
of the quantum state, which are gained through these probes. 

4T. Density ripples. - As we have discussed above, fluctuations play a central role 
in the physics of ID Bose gases. It is thus essential that our method allows the probing 
of a single realization of a ID Bose gas. In this way, repeating the experiment many 
times not only gives access to the dynamics but also to the statistical distribution of the 
fluctuations. It is thus possible to obtain a much deeper insight into the quantum states 
than would be possible if only mean values of observables could be measured. 

A single quasi-condensate that is released and expands in TOF forms strong density 
speckles along the ID axis (see figure 3a). These speckles are a direct consequence 
of the fluctuating phase in the trapped system. In fact, the corresponding gradient 
V6*(z) can be interpreted as a velocity field. In expansion this stochastic velocity field 
transforms into position space realizing a characteristic speckle pattern atop the average 
density profile. Analyzing the correlations in these patterns and comparing them to 
simulated results obtained from an Ornstein-Uhlenbeck stochastic process allows us to 
determine the temperature of the gas [49, 50] as shown in figure 3b. This is a powerful 
tool that works as well for 2D systems [51, 52]. In the experiments it is primarily used 
to characterize the initial gas before the splitting. However, it can also be used for the 
study of the evaporative cooling process [53, 54] or thermalization (see section 6 ). 

4’2. Phase correlation functions. - The interference pattern of two quasi-condensates 
as depicted in the lower panel of figure 2 provides a powerful probe for the dynamics 
of the system. In our case the relative phase fluctuates along the length of the system. 
In general, the position of the fringes in an interference pattern is determined by this 
relative phase between the two interfering waves. The meandering fringe pattern in the 
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Fig. 2. - (a) Overview of the available probes in our setup. The transversal probe is primarily 
used to measure temperature by extracting the density ripple patterns forming in TOF (see 
section 4T). The vertical probe images the full matter-wave interference pattern containing the 
entire relative phase field of the two gases (see section 4’3 and 4'2). The longitudinal probe 
records the interference pattern integrated along the ID direction. It can also be used to measure 
the number balance by intentionally separating the clouds using a momentum transfer during 
the trap switch-off. (b) Examples of interference patterns measured with the vertical imaging 
system right after the splitting {t = 0 ms) and after time evolution {t > 0 ms). The bending of 
the fringes reflects the randomization of the relative phase during the dynamics. Figure adapted 
from [16]. 


images thus directly reflects the local in situ relative phase, which can be reconstructed 
by fitting each local pixel row in the interference pattern with a sinusoidal function. 

Right after the splitting the two halves of the system are almost perfectly phase 
correlated as the shot noise energy is introduced only into the density quadrature, but 
not the phase. The relative phase is almost zero and the fringes are straight. Over time 
this coherence is lost and the fringe patterns become more random. This coherence is 
due to a dephasing of the phononic modes in the relative degrees of freedom. To analyze 
this process it is instructive to study the correlation function of the relative phase field 


( 12 ) 


C{z,z') 
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Fig. 3. - In a TOF measurement the in-situ phase fluctuations transform into density speckle 
patterns (a). The correlations in these patterns are used to extract the temperature by fitting 
them with simulated data (b). The insets show typical density ripple patterns with the displayed 
correlations. Figure adapted from [16]. 


Here, correspond to the field operators of the left and right gas and z and z' are 
two points along the axial direction of the system. In the last step we have assumed that 
density fluctuations can be neglected, which is a very good approximation in the quasi¬ 
condensate regime. In the experiment, the expectation value is realized by averaging 
over many identical realizations. 

For the coherent phase field right after splitting the correlation function is close to one 
over all relative distances z = z — z'. After approximately 15 ms the systems settles into 
a steady state, where correlations decay exponentially with z. For a ID Bose gas this ex¬ 
ponential decay corresponds to thermal correlations, with the characteristic length scale 
of the decay A being directly related to the temperature T via A = h^nid/mkBT. How¬ 
ever, while showing characteristic thermal-like correlations, the relaxed state is markedly 
different from thermal equilibrium, as its temperature ksTeg = gnidj^ can be identified 
with the shot noise energy that was introduced during the splitting process. Is is thus 
significantly smaller than the initial temperature T of the system. At the same time, 
the common degrees of freedom still show a temperature comparable to T. The system 
has thus not fully thermalized, but rather reached a prethermalized state [55, 56], where 
it already exhibits certain thermal-like features like a temperature. The physical reason 
behind this is that common and relative degrees of freedom fully decouple in the low- 
energy limit for a balanced splitting. No energy can be exchanged so that the system 
can never fully forget its initial state. 

Microscopically this dephasing process can be well understood within the Luttinger 
description. All energy is initially stored in the density quadrature and all phonons are 
initialized in phase. During the time evolution the energy of each mode oscillates be¬ 
tween density and phase with the momentum-dependent frequency ujk, which eventually 
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leads to a dephasing. The thermal nature arises from the occupations of the modes. 
Because of the linear dispersion relation we find that the splitting prepares the relative 
degrees of freedom with occupation numbers that decay as 1/fc for increasing momentum 
k. All modes thus obtain the same amount of energy from shot noise, which, after de¬ 
phasing, makes the state indistinguishable from a thermal state with the corresponding 
temperature. 

More insights can be obtained by studying the details of the correlation functions 
during the relaxation process. Their evolution is plotted in figure 4a [57]. For a given 
point in time the correlations decay exponentially up to a certain crossover distance Zc 
beyond which the long-range order of the initial state prevails. The evolution of this 
crossover point plotted in figure 4b is linear, revealing that the exponentially decaying 
correlations spread through the system in a light-cone-like dynamic with a characteristic 
velocity. This process is driven by the dephasing of the phononic modes of the initial 
state. Short wavelength modes dephase faster than long wavelength modes leading the 
characteristic spread. The velocity can be identified with the speed of sound of the 
phonons, which thus act as carriers of information in the system. This observation 
provides a direct connection between the establishment of thermal properties and the 
propagation of correlations in a quantum many-body system. The underlying principles 
are even more general and also govern the distribution of entanglement, with profound 
implications, e.g. for quantum information science and computer simulations of complex 
materials [58, 59, 60]. 

4’3. Full distribution functions. - Another powerful technique to analyze the correla¬ 
tion properties during the relaxation dynamics and especially in the steady state is the 
full distribution function (FDF) of the interference contrast. To introduce the contrast 
as an observable we define the operator 

/•^/2 

(13) A{L)= / dz^{^ {z,t)^r{z,t), 

Jl/2 

which corresponds the interference term of the bosonic field operators integrated over a 
length L. The magnitude of A{L) is related to the integrated contrast of the interference 
patterns {C‘^{L)) = {\A{L)\‘^) . Experimentally the distribution of the squared 
contrast normalized by the mean squared contrast a = C‘^/{\C\^) is less prone to sys¬ 
tematic errors and therefore favorable. Recording the shot-to-shot fluctuations of this 
quantity gives us the full distribution function P{a)da of the probability to observe a 
contrast in the interval a + da. The FDFs therefore contain the information about all 
even moments of the interference operator (13) defined above 

ql|2m\ foo 

(14) J ' = (a™) = / P{a)a^da. 

(|A|2)- 'Wo 

Thus, they contain much more information about the quantum state than the two-point 
correlation function introduced earlier. 
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Fig. 4. - (a) Measured phase correlation functions (circles) of the evolution following the splitting 
together with the Luttinger liquid predictions (solid lines) taking the trap as well as the finite 
optical resolution into account. The color encodes time going from blue (1 ms after the splitting) 
to red (9 ms). The green exponential curve is the prediction for the final prethermalized steady 
state, (b) Evolution of the crossover distance Zc between the exponentially decaying correlations 
and the plateau with long range order. The linear behavior shows that the thermal correlations 
appear locally and spread through the system in a light-cone-like fashion. Figure adapted 
from [16]. 


Figure 5a shows the evolution of the mean squared contrast as a function of time, 
Figure 5b a comparison of the FDFs of the prethermalized state discussed in section 4'2 
and the predictions of the Luttinger liquid model. In figure 5c the FDFs of a system of 
two independent condensates in thermal equilibrium are plotted for comparison. Due to 
the low effective temperature of the prethermalized state its distributions are peaked over 
long integration lengths while the much hotter thermal state in figure 5c exponentially 
decays over all observed length scales. This illustrates the fact that the steady state 
reached after splitting is not the thermal equilibrium of the system. 























12 


B. Rauer, T. Schweigler, T. Langen, J. Schmiedmayer 



evolution time t (ms) 


0 


0 


FULL CLOUD 60 (jm 38 Mtn 27 pm 22 pm 




i 


s: 







01230123012301230123 


Normalised Squared Contrast, 


Fig. 5. - Contrast dynamics and full distribution functions of a coherently split ID Bose gas. 
(a) Measured values of the mean squared contrast for various integration lengths L (points). 
From top to bottom: L = 18, 40, 60,100 pm. The lines show the results of a Luttinger liquid 
calculation for these integration lengths, (b) Full distribution functions after relaxation to 
the prethermalized state. The solid red lines show theoretical equilibrium distributions with 
an effective temperature of Teff = 14: uK, which is significantly lower than the true initial 
temperature of the gas (T = 120 uK). The prethermalized nature of the state is clearly revealed 
by comparing it to the vastly different thermal equilibrium situation shown in (c), which can 
be prepared by creating two completely independent ID Bose gases. Figure adapted from 
Refs. [55, 56]. 


5. Generalized Gibbs ensemble 

The fact that the phonon occupations of the system are preserved during the dynamics 
is deeply rooted in the integrability of the underlying model. Each relative mode acts 
like a harmonic oscillator that does not interact but dephases with respect to the rest of 
the system. This is a general feature of an integrable quantum system, where multiple 
non-trivial quantities are conserved, severely restricting the system’s dynamics. This was 
strikingly visualized in a landmark experiment by Kinoshita et al. [26], which realized 
the quantum analog of the well known (classical) Newton cradle. Even after thousands 
of collisions between its constituents such a system will not reach thermal equilibrium, 
simply because the momenta are conserved and can thus never reach the values given by 
the Bose-Einstein distribution. 

Nevertheless, it has been conjectured that such systems still relax to a maximum en¬ 
tropy state which is given by the density matrix of a so-called generalized Gibbs ensemble 
(GGE) [37] 

(15) 

Here, Z is the partition function, Ij are the operators of the conserved quantities and Xj 
the corresponding Lagrange multipliers. If only energy is conserved this density matrix 
reduces to the well-known canonical or Gibbs ensemble, with temperature being the only 
Lagrange multiplier. If many more conserved quantities exist like the phonon occupations 
in the Luttinger liquid model, many generalized temperatures, one for each conserved 
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quantity are necessary to maximize entropy. 

In our case the occupation numbers of all modes are the conserved quantities. How¬ 
ever, the prethermalized state that we have studied so far is a special case of this ensemble, 
as all temperatures are identical due to the equipartition of energy during the splitting 
process. To demonstrate the presence of a GGE it is thus necessary to change the splitting 
process, such that different modes exhibit different temperatures. The results are shown 
in Fig. 6. Again, the relative phase correlation function can be used to characterize the 
dynamical states of the system. While we were previously showing only one coordinate 
of this function, plotting the full function provides straightforward insights into the new 
occupation numbers. The correlation functions show a trivial maximum on the diagonal 
(zi = Z 2 ) which arises due to the fact that every point is perfectly correlated with itself. 
However, a second maximum arises on the anti-diagonal (zi = —Z 2 ), indicating that 
points that are located symmetrically around the center of the system are more strongly 
correlated. In a simplified model, this implies that modes which are symmetric around 
the center are more strongly occupied than modes with are anti-symmetric around the 
center. A more detailed analysis of the relaxed state allows to extract all mode occu¬ 
pations that are necessary to describe the state [61]. Given these extracted occupation 
numbers the dephasing model also provides a detailed description of the dynamics, which 
proves that the conserved quantities were indeed set during the splitting process. 

Most importantly, these observations visualize, both experimentally and theoretically, 
how the unitary evolution of our quantum many-body system connects to a steady state 
that can be described by a thermodynamical ensemble. 



Fig. 6. - Relaxation dynamics of a coherently split ID Bose gas with different populations for 
different modes. Two-point correlation functions C(z, z') for increasing evolution time, showing 
maxima on the diagonal and the anti-diagonal. The experimental observations (top row) are in 
very good agreement with the theoretical model (bottom row) demonstrating the presence of 
many different temperatures in the system. Figure adapted from [16, 61]. 







































14 


B. Rauer, T. Schweigler, T. Langen, J. Schmiedmayer 


6. Dynamics beyond prethermalization 

In sections 4 and 5 of these notes we demonstrated that the ID Bose gases realized in 
experiment do not relax to thermal equilibrium but to a prethermalized state that can 
be described by a generalized Gibbs ensemble. This behavior roots in the integrability 
of the Lieb-Liniger model and its low-energy approximation, the Luttinger liquid model. 
However, the ID Bose gas realized in our experiments is only nearly-integrable. On the 
one hand radial excited states can affect the ID dynamics and on the other hand the 
harmonic trap breaks integrability of the Lieb-Liniger model (while integrability is still 
retained in the trapped Luttinger liquid model [62]). 

It has been conjectured that in this case the observed prethermalized state is only an 
intermediate steady state on the way to thermal equilibrium, its lifetime being directly 
related to the degree of integrability breaking [63, 64]. The analysis of this scenario in 
the context of classical mechanics has culminated in the important Kolmogorov-Arnold- 
Moser (KAM) theorem [65]. No complete analogue of this theorem has so far been 
found in quantum mechanics [66]. Alternatively, also other behaviour has been sug¬ 
gested, namely that the quasi-particles of the experimentally realized ID Bose gas could 
be unaffected by the radial states [67], leaving the gas fully integrable. Experimental 
investigations into this effect are ongoing in our and other groups [68]. 

However, even within the coherent dynamics the long-term evolution of the system is 
expected to show a rich variety of effects, which we will discuss in the following. 

6T. Recurrences. - We have shown in the previous chapters that the unitary quantum 
evolution of a ID Bose gas can lead to the establishment of thermal properties. This 
does not mean that a true thermal state was reached, but rather that the expectation 
values of certain observables became indistinguishable from the corresponding thermal 
values. In this way the predictions of statistical and quantum mechanics are reconciled. 

However, in a finite system as the trapped system we are dealing with, the unitarity is 
still expected to result in observable consequences as it forces the dynamics to be periodic. 
The important question is how long the timescale of this periodic behaviour will be. In 
the context of our experiments periodic behavior would correspond to a rephasing of the 
phonons (and thus a reestablishment of coherence) after a finite time, which would be 
observable as a phase correlation function close to one, C{z) = 1, over all distances z. 

In a homogeneous system the time between these recurrences can be estimate as Gee = 
L/2c, which corresponds to twice the time to reach the perfectly dephased prethermalized 
state. For typical parameters Gee ~ 30 ms. Surprisingly no signs of these recurrences are 
observed in experiment. The reason for this lies in the mode structure of the trapped 
system. While in the homogeneous case the mode energies are equally spaced Wfe = cfc, the 
modes in a harmonically trapped condensate are described by Legendre polynomials [69]. 
This leads to the modified dispersion relation 

(16) UJj = WaxVj(j + l)/2, 

where Wax is the trap frequency of the axial harmonic confinement and j is the mode 
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index. For the given parameters the incommensurate mode frequencies shift the first 
significant revival in the trapped case to about 200 ms, which is challenging to study 
in experiment [62]. Figure 7 shows a comparison of the phase correlation dynamics 
after splitting for the homogeneous and the trapped case. While the initial dephasing 
dynamics is very similar in both traps the revival structure is quite different, as expected 
from the dispersion relations. A classical analogy for these dynamics is the behavior of a 
collection of uncoupled pendula of different length, which only rephase if their frequencies 
are commensurate. 

To actually measure recurrences in experiment it would be beneficial to trap atoms 
in a box shaped potential. Flat bottom traps were recently realized for 3D and 2D 
systems [70, 71]. They are a powerful tool to investigate Bose-Einstein condensation or 
the Kibble-Zurek mechanism unperturbed by trap effects. 

6’2. Imbalanced splitting. - Another relaxation mechanism that is captured by the 
low-energy description is the dephasing due to imbalances in the splitting process. In 
practice, the two wells of the double-well potential can never be perfectly balanced during 
the splitting process. This leads to relative fluctuations of the overall number of atoms in 
each well. The gas which ends up with more atoms is characterized by a slightly higher 
chemical potential and speed of sound. These relative differences lead to a dephasing 
of the two gases with a characteristic velocity c' = (c/ — Cr)/2 [16, 47]. If the atom 
number difference between the two gases is very small this process will thus be much 
slower than the initial relaxation to the prethermalized state. However, on long time- 
scales it will lead to a state in which common and relative degrees of freedom share the 
same temperature. For an observer, the state will thus be indistinguishable from thermal 
equilibrium, highlighting again the importance of dephasing and the role of observation 
for the understanding of thermalization. 

7. Application: Interferometry with squeezed states 

We will end with an application of the well-characterized matter-wave interferometer 
that we have introduced during the course of these notes. With this we aim to indicate 
how the fundamental research on non-equilibrium dynamics might also have immediate 
technological impact in the near future. 

The binomial splitting of a single gas that we discussed above is only a good ap¬ 
proximation in the limit of non-interacting atoms. Even for weak interactions, as the 
ones in the present case for ®^Rb, the splitting has to be very fast to reach the binomial 
splitting limit. For a slower splitting, interactions will start to play a role and lead to 
the development of correlations between the atoms. These correlations are a valuable 
resource for precision measurement devices. 

Experimentally, the speed of the splitting can easily be controlled using the atom 
chip. While an infinitely fast splitting leads to a relative atom number variance of N, 
the variance resulting from a splitting taking a finite time is reduced by the so-called 
squeezing factor , where ctjv is the standard deviation of the relative number 
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separation z (micron) 



Fig. 7. - Time evolution of the relative phase correlation function for the homogeneons (left, a 
and c) and trapped (right, b and d) systems. The color-scale indicates the degree of correlation 
(red: high correlation, blue: low correlation). The top row illustrates the relaxation to the 
prethermalized state. In the homogeneous case, the initial state is re-established at times which 
are multiples of the system length divided by the characteristic velocity. In the trapped case, 
the recurrences are only partial and the more complex structure is due to the incommensurate 
ratios of the mode frequencies. In this time window (0 — 300 ms), the strongest recurrence is 
observed at 202ms (Here, tUax = Ztt • 7Hz). Reproduced from [62] 


distribution. The slower the splitting the lower the factor and thus the stronger the 
squeezing. The corresponding spin squeezing factor /{cos(p) {4> again denoting 

the relative phase) can be understood as an entanglement witness, i.e. an observable 
that signals the presence of genuine multi-particle entanglement [72]. The presence of 
this entanglement in the states created by the splitting leads to a gain in measurement 
precision, which can not be achieved with classical states [73]. 
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To actually utilize number squeezing in our setup we need to devise a full interferomet¬ 
ric sequence. Apart form the splitting process already described at length we need two 
further building blocks to achieve this. The first is a mechanism to introduce a relative 
phase shift between the two arms of the interferometer, which emulates the measurement 
signal in a possible application. Experimentally we realize this by tilting the double well 
potential after splitting so that the two gases experience a different gravitational poten¬ 
tial and accumulate a phase difference. Varying the time the system is kept in this 
state controls the overall phase shift. As a second building block, we need to employ a 
recombiner that allows for measurements of the relative atom number in the two arms of 
the interferometer (in analogy to the second beam splitter in an optical Mach-Zehnder 
interferometer). This can be achieved by accelerating the two gases onto each other while 
keeping a barrier between them that is small enough to allow for inter-well tunneling. 
In this process the relative phase is mapped to a relative population difference, just like 
in the case of two wave packets that simultaneously impinge on a semi-reflective barrier 
from different sides. In addition, the relative phase (j) between the two wells can also 
be measured using the standard matter-wave interference procedure that was already 
described in section 4. 

Figure 8 shows the experimentally observed population imbalance of a squeezed initial 
state as a function of the phase accumulation time The fringe contrast of the average 
values is damped due to phase diffusion. Naively, this phase diffusion would be expected 
to be much more severe, but the presence of a long-lived prethermalized state limits 
its deteriorating effects. Comparing the observed decay time to the one expected for a 
coherent state of = 1 (dashed line) illustrates the gain in interferometric precision 
when using a squeezed input state. The best spin squeezing achieved in this setup is 

= —7.8 ± 0.8 dB [74], corresponding to genuine multi-particle entanglement of 150 
atoms. This result could in the future be increased by the use of optimized splitting 
ramps [75], and outlines the way for interferometric sensing of local forces in atom chip 
configurations. 

8. Conclusion 

The relaxation of isolated quantum many-body systems is a major unsolved problem 
connecting statistical and quantum physics. Understanding such relaxation processes 
remains a challenge despite considerable efforts. Experiments with ultracold quantum 
gases (in general) and ID Bose gases (in particular) allow the realization and manipu¬ 
lation of well-controlled and truly isolated quantum systems. As we have shown, this 
provides unique opportunities to study and understand non-equilibrium phenomena. For 
example, the results discussed in these notes demonstrate for the first time several char¬ 
acteristic aspects of these dynamics, including the existence of a stable, thermal-like 
prethermalized state and its dynamical, light-cone-like emergence. Furthermore, the 
connection of the prethermalized state with generalized statistical ensembles, and thus 
of the unitary quantum evolution and statistical mechanics was highlighted. The progress 
in this field is rapid and we expect it to continue to have profound implications for our 
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Phase accumulation time (ms) 

Fig. 8. - Output signal of the integrated Mach-Zehnder Interferometer. The normalized popula¬ 
tion difference z = Ni — Nr/{Ni -|- Nr) = n/Nt between the two wells is measured as a function 
of time It exhibits interference fringes and a damping due to phase diffusion. Grey dots: 
imbalance of individual experimental realizations; black dots: ensemble average ( 2 ); red curve: 
theoretical prediction taking into account phase diffusion; dashed black line: expected signal for 
a classical coherent state without squeezing. Reproduced with permission from [74]. 


understanding of isolated quantum many-body systems. 
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